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Abstract 



Inclusive photoproduction of D*^ in ep collisions at HERA has been measured 
with the ZEUS detector for photon-proton centre of mass energies in the range 
115 < W < 280 GeV and photon virtuality Q 2 < 4GeV 2 . The cross section 
a ep -+D*x integrated over the kinematic region p®* > 3 GeV and — 1.5 < rj D * < 1.0 
is (10.6 ± 1.7(stat.) zbjj (syst.)) nb. Differential cross sections as functions of p^*, 
r] D and W are given. The data are compared with two next-to- leading order per- 
turbative QCD predictions. For a calculation using a massive charm scheme the 
predicted cross sections are smaller than the measured ones. A recent calculation 
using a massless charm scheme is in agreement with the data. 
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1 Introduction 



Production of the heavy quarks c and b at HERA is dominated by photoproduction, where 
a quasi-real photon with negative four momentum squared, Q 2 , close to zero is emitted 
by the incoming electron and interacts with the proton. Heavy quark photoproduction 
can be used to probe perturbative QCD calculations with a hard scale stemming from 
the heavy quark mass and the high transverse momentum of the produced parton. At 
leading order (LO) in QCD two types of processes are responsible for the production of 
heavy quarks: the direct photon processes, where the photon participates as a point-like 
particle which interacts with a parton from the incoming proton, and the resolved photon 
processes, where the photon is a source of partons, one of which scatters off a parton from 
the proton. Charm quarks present in the parton distributions of the photon, as well as 
of the proton, lead to processes like eg —>■ eg, which are called charm flavour excitation. 
In next-to-leading order (NLO) QCD only the sum of direct and resolved processes is 
unambiguously defined. Two types of NLO calculations using different approaches are 
available for comparison with measurements of charm photoproduction at HERA. The 
massive charm approach ^ |3[] assumes light quarks to be the only active flavours 
within the structure functions of the proton and the photon, while the massless charm 
approach M 0] also treats charm as an active flavour. 

The total charm photoproduction cross section a 1P ^ C cx has recently been measured 
at HERA H at a -yp centre of mass energy W ~ 200 GeV, and was found to be about 
one order of magnitude larger than fixed target data. These measurements were compared 
with NLO QCD calculations and with a calculation based on a semihard approach of 
QCD ||. The precision of the comparison was limited by large systematic uncertainties 
in both data and theory. In the case of the data the uncertainty was due to the necessity 
of extrapolating the measured cross sections into kinematic regions not accessible to the 
experiments. The uncertainties of the theoretical predictions were generated by the de- 
pendence of the NLO calculations on the charm mass, m c , the QCD renormalization and 
factorization scales, and on the parton density parametrizations assumed for the proton 
and the photon. 

In this study we use a sample of D*(2010) ± mesons collected with the ZEUS detector 
during 1994. The sixfold increase of the data sample with respect to our previous analy- 
sis || allowed the measurement of the differential distributions in p^ , i] D and W, where 
p±* is the D* transverse momentum with respect to the beam axis and r] D * is the pseu- 
dorapidityQ of the D*. The measurement of the ep — > D*X cross sections was performed 
in the restricted kinematic region Q 2 < 4 GeV 2 , 115 < W < 280 GeV, p®* > 3 GeV and 
— 1.5 < i] D * < 1.0. The results are compared with NLO QCD predictions calculated in 
both the massive and the massless charm approach. The theoretical uncertainties in the 
calculations are strongly reduced in the restricted kinematic region, thus allowing a more 
precise comparison of the perturbative QCD calculations with our data. The results of 

1 The pseudorapidity rj is defined as — ln(tan |), where the polar angle 9 is taken with respect to the 
proton beam direction. 
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the calculations using a massless charm scheme are sensitive to the charm content of the 
photon and are insensitive to that of the proton [f|. [| . 

D* mesons are reconstructed from their decay products through the two decay modes^: 

D* + - D\+ -> (K-n + )n+, (1) 
D *+ ^ D \+ ^ ( j ^- 7r + 7r + 7r -) 7r +. ( 2 ) 

The small mass difference M(D*) - M(D°) = 145.42 ± 0.05 McV [|Kj yields a low mo- 
mentum pion ("soft pion", its) from the D* decay and prominent signals just above the 
threshold of the M(Knirs) — M(Kir) and M^Kim-Kits) — M(Kir7nr) distributions, where 
the phase space contribution is highly suppressed |TT . 



2 Experimental Conditions 

The data presented in this analysis were collected during the 1994 running period using 
the ZEUS detector at HERA, where a positron beam with energy E e = 27.5 GeV collided 
with a proton beam with energy E p = 820 GeV. A total of 153 colliding bunches were 
stored in HERA, together with additional 17 proton and 15 positron unpaired bunches 
intended for studies of beam-induced backgrounds. The time interval between bunch 
crossings was 96 ns, and the typical instantaneous luminosity was 1.5T 30 cm~ 2 s _1 . The 
r.m.s. of the vertex position distribution along the beam direction was 12 cm. The total 
integrated luminosity used in this analysis is 2.99 pb -1 . 



2.1 The ZEUS Detector 

A detailed description of the ZEUS detector can be found in ref. |12|, [L3|]. Here we present 
a brief description of the components relevant to the present analysis. 

Charged particles are measured by the Central Tracking Detector (CTD) [|14[] which 
operates in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD is a drift chamber consisting of 72 cylindrical layers, arranged in 9 superlayers. 
Superlayers with wires parallel to the beam axis alternate with those inclined at a small 
angle to give a stereo view. The single hit efficiency of the CTD is greater than 95% and 
the measured resolution in transverse momentum for tracks with hits in all the superlayers 
is a p± /p ± = 0.005p ± ©0.016 (p ± in GeV). 

Surrounding the solenoid is the uranium-scintillator calorimeter (CAL) Ul5| , which 
is divided into three parts: forward,]] barrel and rear covering the polar regions 2.6° to 

2 In this analysis the charge conjugated processes are also included. 

3 Throughout this paper we use the standard ZEUS right-handed coordinate system, in which X — 
Y = Z = is the nominal interaction point, the positive if-axis points in the direction of the protons 
(referred to as the forward direction) and the X-axis is horizontal, pointing towards the centre of HERA. 
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36.7°, 36.7° to 129.1° and 129.1° to 176.2°, respectively. The CAL covers 99.7% of the 
solid angle, with holes of 20 x 20 cm 2 in the centres of the forward and rear calorimeters 
to accommodate the HERA beam pipe. Each of the calorimeter parts is subdivided 
into towers which are segmented longitudinally into electromagnetic (EMC) and hadronic 
(HAC) sections. These sections are further subdivided into cells each of which is read out 
by two photomultipliers. From test beam data, energy resolutions of <Je/E = 0.18 /VE 
for electrons and ge/E = 0.35/v^ for hadrons (E in GeV) have been obtained. The 
timing resolution of a calorimeter cell is less than 1 ns for energy deposits greater than 
4.5 GeV. In order to minimise the effects of noise due to the uranium radioactivity on the 
measurements, all EMC (HAC) cells with an energy deposit of less than 60 (110) MeV are 
discarded from the analysis. For cells without energy deposits in neighbouring cells this 
cut was increased to 80 (120) MeV. 

Proton-gas events occurring in front of the nominal ep interaction region are out of 
time with respect to the ep interactions and may thus be rejected by timing measurements 
made by the CAL and by scintillation veto counter arrays located at Z — —730 cm, 
Z = —315 cm and Z = —150 cm. 

The luminosity was determined from the rate of the Bethe-Heitler process e + p — > e + ~/p, 



where the photon is measured by a calorimeter |L6[ located at Z — —107 m in the HERA 
tunnel in the direction of the positron beam. 



2.2 Trigger 

The ZEUS detector uses a three level trigger system |T2J. For the analysis presented 
in this paper the following trigger strategy was chosen to identify D* candidates in the 
central region of the detector. 

In the first level trigger the calorimeter cells were combined to define regional and 
global sums which were required to exceed one of the following thresholds: 

• Total energy greater than 15 GeV, 

• Total EMC energy greater than 10 GeV, 

• Total transverse energy greater than 11 GeV, 

• EMC energy in the barrel calorimeter greater than 3.5 GeV, 

• EMC energy in the rear calorimeter greater than 2 GeV. 



In addition, at least one CTD track coming from the ep interaction region was required. 
Events with timing measured by the veto counters, inconsistent with an ep interaction, 
were removed. 
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In the second level trigger beam-gas events were rejected by exploiting the excellent 
timing resolution of the calorimeter. Also events were rejected in which the vertex deter- 
mined by the CTD was not compatible with the nominal ep interaction region. 

In the third level trigger (TLT) the full event information was available. Calorime- 
ter timing cuts were tightened for further rejection of the remaining beam-gas events. 
Events were required to have a transverse energy outside a cone of 9 = 10° with respect 
to the proton direction, E^ 10 , above 12 GeV. The online value overestimates the offline 
£j0>io ^ ue £ G simplified event reconstruction at the TLT. In addition, the following 
requirements were made: 

• p™ ax+ > 0.5 GeV, 

• p™ ax ~ > 0.5 GeV and 

• p™ ax+ + p™ ax ~ > 2.0 GeV, 



where the quantity p™ ax+ (p™ ax ~) is defined as the maximum transverse momentum of 
any positive (negative) track associated with the reconstructed event vertex in the polar 
angle range 15° < 6 < 165°. 



3 Data Analysis 

3.1 Offline Data Selection 

The event sample accepted by the trigger algorithm was processed using the standard 
offline ZEUS detector calibration and event reconstruction code. To define an inclusive 
photoproduction sample, the following requirements were imposed: 



A reconstructed vertex with at least three tracks associated to it. 
E{> 10 ° > 12 GeV. 

No scattered positron found in the CAL according to the algorithm described in 
ref. |T7[]. This requirement removes deep inelastic scattering (DIS) neutral current 



events, thereby restricting Q to below 4 GeV . The corresponding median Q is 
about 5-10" 4 GeV 2 . 



0.1 < yjB < 0.7. Here yjs is the Jacquet-Blondel |18| estimate of y, the fraction of 
the positron energy carried by the photon in the proton rest frame. It is defined as: 



y^B 



T,j(E - p z 
2E P 
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where the sum runs over all calorimeter cells and p z is the Z component of the 
momentum vector assigned to each cell of energy E. The lower y^ cut rejects events 
from a region where the acceptance is small because of the trigger requirements. 
The upper cut rejects possible background from DIS events in which the scattered 
positron has not been identified, and therefore is included in the y^ calculation, 
thus producing a value of y^B closer to 1. 

With these requirements, an inclusive photoproduction sample of about 450,000 events 
was selected. The jp centre of mass energy of these events was calculated from ?/jb via 
the expression Wjb = ^ Ay^E p E e and ranges from 100 to 250 GeV. A systematic shift 
observed in the reconstructed values of Wjb with respect to the true W of the event, due 
to energy losses in inactive material in front of the calorimeter and particles lost in the 
rear beam pipe, was corrected using Monte Carlo (MC) techniques || [17|]. The centre of 
mass energy range covered by the photoproduction sample is then 115 < W < 280 GeV. 

A D* reconstruction algorithm was applied to all selected events. This algorithm 
combines the reconstructed tracks in each event to form D* candidates assuming the 
decay channels (HD or It uses the mass difference technique to suppress the high 
background due to random combinations from non-cc events, which have a much higher 
cross section. Only tracks associated with the event vertex and having p± > 0.2 GeV and 
\r]\ < 1.75 are included in the combinations. 

3.2 The D* + — > D°iz^ — > (K-tz + )tz^ Decay Channel 

In each event tracks with opposite charges and p± > 0.5 GeV were combined into pairs 
to form D° candidates. The invariant mass M(Ktt) of each pair was calculated. No 
particle identification was used, so kaon and pion masses were assumed in turn for each 
particle in the pair. A third track, assumed to be the soft pion, its, with a charge 
opposite to that of the particle taken as a kaon, was then added to the combination. The 
mass difference AM = M(K7nrs) — M(Ktt) was evaluated. Only Ktttt,s combinations 
with —1.5 < f] < 1.0, for which the acceptance is high, were kept. As a result of the 
£t<9>io > Yi GeV cut the acceptance of the KniTs combinations having p± below 3 GeV 
is very small. Thus pj_(KinTs) > 3 GeV was required. This cut excludes almost all soft 
pions with p±(7Ts) < 0.2 GeV, so the overall cut of p± > 0.2 GeV on all tracks causes 
essentially no loss of D* candidates. 

In Fig. la the AM distribution for combinations with 1.80 < M(Ktt) < 1.92 GeV 
is shown. A clear peak at the nominal value of M(D*) — M(D°) [|10J is observed. 
To determine the combinatorial background under the peak, combinations in the range 
1.68 < M(Ktt) < 2.04 GeV were used, in which both tracks forming the D° candidates 
had the same charge (these will be referred to as wrong charge combinations). The AM 
distribution from these combinations, normalized to the number of right charge combina- 
tions in the range 155 < AM < 180 MeV, is also shown in Fig. la. 
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The M(Ktt) spectrum corresponding to combinations having a mass difference in the 
range 143 < AM < 148 MeV is shown in Fig. lb. A clear peak at the nominal D° 
mass of 1.8645 ± 0.0005 GeV jUj is observed. The combinatorial background was again 
determined by using wrong charge combinations. The M(Ktt) distribution from these 
combinations was normalized to the number of right charge combinations in the range 
2.0 < M(Ktt) < 2.5 GeV. A clear excess of right charge combinations with respect to 
the combinatorial background was observed for masses below that of the D° meson. This 
excess is reproduced by the MC simulation (when a sample including all decay modes of 
the D° is used), and is mainly due to events in which a D° decays into Kinr® and the 
extra 7r° is not included in the D° invariant mass reconstruction. 

The signals from Fig. 1 were fitted with the maximum likelihood method to a sum 
of a Gaussian (describing the signal) and a functional form (describing the background 
shape) of A ■ (AM — m n ) B for Fig. la and exp(C + D ■ M(Ktt)) for Fig. lb. The mass 
and width values obtained were: AM = 145.44 ± 0.08 MeV, a = 0.68 ± 0.08 MeV from 
Fig. la and M(D°) = 1860 ± 4 MeV, a = 29 ± 4 MeV from Fig. lb. The mass values 



obtained are in agreement with the PDG values fljj| . The width values agree with our 
MC simulations. The observed signals confirm that the reconstructed D* candidates 
come from channel ([!]). The contribution of other D° decay modes to this signal is 
negligible as estimated by MC studies. The number of reconstructed _D*'s was determined 
by subtracting the normalized background distributions described above from the right 
charge distributions. The use of a wider M(Kir) range for the wrong charge combinations 
compared with the signal region reduces the statistical error of the subtraction procedure. 
After subtracting the background from the AM distribution of Fig. la, a signal of 152 ±16 
reconstructed _D*'s was obtained, consistent with the corresponding subtracted signal from 
the M(Ktc) distribution. 



3.3 The D* + — > D°7Zg — > (K-tt+tz+tz-)^ Decay Channel 

For this channel, four tracks were combined to form a D° candidate. The invariant 
mass M(Ktt7ttt) was calculated for combinations having a total charge of zero. Since 
no particle identification was used, kaon and pion masses were assumed in turn for each 
particle in the combination. A fifth track, assumed to be the soft pion, tc,s, with a 
charge opposite to that of the kaon was added to the combination and the mass difference 
AM = M(Kn7nT7Ts) — MIKtttttt) was determined. Minimum transverse momenta of 
0.5 GeV for the track taken to be the kaon and of 0.2 GeV for the other tracks were 
required. 

Due to the higher number of decay particles in this channel, the combinatorial back- 
ground is higher than that of channel ([[]). In order to achieve an improved signal to 
background ratio, we required the transverse momentum of the K r K r n / n r Ks combination to 
be above 4 GeV, and to satisfy the condition p^- 7T7Tnn s > g 2 .^ >10 °. This latter require- 
ment removed about one third of the combinatorial background, i.e. the contribution 
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from events with high E 9 ^ 10 " in which a combination with relatively low transverse mo- 
mentum was found. No reduction in the number of signal events was observed. Finally, 
the reconstructed D* candidates were required to be in the same pseudorapidity range as 
for channel ([!]), i.e. —1.5 < i] < 1.0. 

Fig. 2a shows the AM distribution for those combinations with 1.81 < M(Kii7ni) < 
1.91 GeV. The M(Kinnr) spectrum for combinations in the range 143 < AM < 148 MeV 
is shown in Fig. 2b. A smaller window around the nominal D° mass was used in this 
channel because of the better mass resolution due to the lower average transverse mo- 
menta of the D° decay particles. The signals from Fig. 2 were fitted to a sum of a Gaus- 
sian (describing the signal) and a functional form (describing the background shape) 
of A ■ (AM -m^) B for Fig. 2a and P 2 (x) = C + Dx + Ex 2 , where x = M(Kititit) 
for Fig. 2b. The mass and width values obtained were: AM = 145.42 ± 0.11 MeV, 
a = 0.83 ± 0.11 MeV from Fig. 2a and M(D°) = 1859 ± 3 MeV, a = 19 ± 3 MeV from 



Fig. 2b. The mass values obtained are in agreement with the PDG values p0| , and the 
widths agree with the MC values. 

There are two main sources of background to the D* signals in this decay channel. The 
first one is the combinatorial background coming from events or tracks in which no D* 
decaying through this channel is produced. To determine this background, wrong charge 
combinations with total charge ±2 for the D° candidate and total charge ±1 for the D* 
candidate were used. The distributions of the wrong charge combinations were normalized 
to the number of right charge combinations for AM in the range 155 < AM < 180 MeV, 
and for M(K7nnr) in the range 2.0 < M(K7nnr) < 2.5 GeV and are shown in Fig. 2. 
After the background subtraction, 199 ± 29 reconstructed D* J s were observed in the AM 
distribution. The second source of background is due to events in which a D* decaying in 
this channel was produced, but more than one combination was reconstructed inside the 
signal region, due to the erroneous assignment of the kaon mass to a pion with the same 
charge from the D° decay. The fraction of D*'s with a wrong (K, it) mass assignment 
was found from MC calculations to be 16%. After subtracting this contribution, a final 
signal of 167 ± 25 reconstructed Z}*'s was obtained. Monte Carlo studies show that the 
contribution of other D° decay modes to this signal is about 3% and was neglected. 



Monte Carlo Simulation 



con- 



The Monte Carlo programs PYTHIA |19| and a recent version of HERWIG |20, |21 
taining an improved description of transverse energy distributions were used to model the 
hadronic final states in cc production and to study the efficiency of the data selection cuts. 
Both programs are general purpose event generators which include QCD LO matrix ele- 
ments for charm photoproduction. They simulate higher order QCD radiation by parton 
shower evolution in the initial and final states, taking into account coherence effects from 
the interference of soft gluon amplitudes. Fragmentation into hadrons is simulated with 



a cluster algorithm [E2J in the case of HERWIG, and with the LUND string model [23 
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in the case of PYTHIA. In HERWIG, the lepton-photon vertex is calculated exactly for 



the direct photon processes and the equivalent photon approximation |23] is used for the 
resolved processes. In PYTHIA the Weizsacker- Williams approximation [25) is used in 
both cases. 

Large samples of cc events were generated with both MC programs. Direct and re- 
solved photon events, including charm excitation, were generated using several parton 
distribution parametrizations for both the proton and the photon. The MRSG [EH] 



parametrization for the proton and the GRV-G HO [^7| f° r the photon were used to 
produce the reference samples. These samples had at least 20 times the statistics of the 
data, so their contribution to the statistical error was neglected. It was found from MC 
studies that, in the kinematic range used, the results are insensitive to contributions from 
charm excitation in the proton. However they are sensitive to charm excitation in the 
photon. The differences between results obtained with and without charm excitation in 
the photon are included in the systematic errors. A large sample of bb events was also 
generated with PYTHIA to allow an estimation of the fraction of photoproduced D*'s 
originating from this process. 

Events containing at least one charged D* decaying into channel ([!]) or (|2|) were pro- 
cessed through the standard ZEUS detector and trigger simulation programs and through 
the event reconstruction package. 



5 Cross Section Determination 



The integrated D* electroproduction cross section in the kinematic region defined by the 
selection cuts Q 2 < 4GeV 2 , 115 < W < 280 GeV, -1.5 < r/ D * < 1.0 and > 3 
or 4 GeV is calculated using the formula: 



N D * 

corr 

L ■ Br 



where N® orr is the acceptance-corrected number of D*, Br is the combined branching ratio 
of a given channel (see Table 1) and C = 2.99 ± 0.05 pb _1 is the integrated luminosity. 

The usual method to correct data by the ratio of generated to reconstructed MC 
events is valid when the MC well describes the data distributions of quantities used 
in the analysis. In our case most of these distributions are well described by the MC 
events. However, for the E e ± 10 distribution the data yield higher values compared to 
the events simulated with PYTHIA. Therefore the reference MC used to calculate the 
acceptance for channel (1) was HERWIG. In the case of channel (0) PYTHIA was used 
since HERWIG does not reproduce the decay width of the resonances, and hence does 
not describe correctly the D° decays into four particles. 

In order to minimize the MC dependence of the corrections for the E e > 10 ° cut in both 
channels, a two step procedure was used to calculate iV^* r : 
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(a) A weighting factor u\ was obtained from the MC simulation and applied to each 
reconstructed D* candidate. These events were corrected for the tracking efficiency 
of the D* reconstruction algorithm and for all trigger and event selection cuts except 
for the i^ >10 ° cut. The factor ui\ is defined as the number of generated D*'s divided 
by the number of reconstructed D*'s. It is calculated in a three-dimensional grid 
in the reconstructed quantities p®* , rf* and Wjb- The average values of uj\ are 
2.99 ± 0.03 and 3.07 ± 0.04 for channel (1) and channel (2), respectively. The 
variation of lo% as a function of any of the three grid variables is less than a factor 
of two in the restricted kinematic range. 

(b) To calculate the correction for the £^ >10 ° > 12GeV cut we used an independent 
D* data sample selected for channel (1) with no cut on the energy deposition in the 
calorimeter at the TLT. Only tracking information was used for this sample selection 
at the TLT. The corresponding correction factor uj 2 is given by the total number of 
_D*'s (corrected as in (a) above) obtained with the TLT tracking selection divided 
by the number of those _D*'s (corrected as in (a) above) in events which satisfy the 
cut E^ 10 > 12GeV. The correction by this method increases the statistical error 
of the result, but reduces the systematic error due to the different MC modeling of 
the E^ 10 distribution. The factor u 2 was applied as an overall weight for the total 
cross section calculation or bin-by-bin in each differential cross section distribution. 
Since the u>2 correction is independent of the tracking acceptance (which is included 
in u>i), the u 2 weights of channel (1) were used also for the N®* rr evaluation of 
channel (2). The average values of u>2 are 1.8 ± 0.2 and 1.4 ± 0.1 for p± > 3 GeV 
and p± > 4 GeV, respectively. 

The reconstructed number of D* mesons after background subtraction (iVj^* as ), the 
acceptance corrected number of produced D*'s [N^ rr = a^u^iV^*^) , the branching ratios 
Br [|K| and the integrated cross sections for both decay channels are presented in Table 1. 
For comparison with channel (2), the value obtained for channel ([I]) requiring p®* > 4 GeV 
is also shown. The cross sections of both channels are in good agreement for the same 
kinematic region. 

The differential cross sections da /dp®*, d<j/dr] D * and da/dW were calculated using 
the same correction procedure in five bins in p®* (3-4; 4-5; 5-6; 6-8; 8-12 GeV), three 
bins in r] D * ((-1.5)-(-0.6); (-0.6)-(0.1); (0.1-1.0)) and three bins in W (115-170; 170-225; 
225-280 GeV). For each variable (p±* , 1] D , W), the other two variables are integrated over 
their kinematic region. The combinatorial background was subtracted bin-by-bin from 
each distribution using the method described in section 3. The dajdp 1 ^ distributions for 
both channels (1) and (2) are shown in Fig. 3. The dajdrf* and da/dW distributions 
for channel (1) are shown in Fig. 4. The results are in agreement with those of the HI 
untagged photoproduction data [0], also shown in Fig. 3, with the same D* rapidity range 
and with 95 < W < 268 GeV, which is quite similar to the kinematic range of this 
measurement. 
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Channel 


(GeV) 


N D * 

meas 


N D * 

corr 


Br 


Oe V ^D*x (nb) 


® 


> 3 

> 4 

> 4 


152 ± 16 
97 ±11 
167 ±25 


828 ± 129 
348 ± 53 
739 ± 128 


0.0262 ±0.0010 
0.0262 ±0.0010 
0.051 ±0.003 


10.6 ± 1.7±i;jj 
4.5 ±0.7 ±0.6 
4.8±0.8±S;g 



Table 1: 

GeV and — 1.5 < r\ 



Cross section a ep ^o*x and related quantities for Q 2 < 4 GeV 2 , 115 < IF < 280 
D * < 1.0 . 



5.1 Systematic Uncertainties 



A detailed study of possible sources of systematic uncertainties was carried out, which we 
now summarize. 

Trigger acceptance uncertainties are mainly due to the different energy distributions 
predicted by the two event generators PYTHIA and HERWIG with and without flavour 
excitation. The uncertainties were determined from the difference between the cross 
sections obtained with both MC generators. For channel (1) the largest shift was —11%. 
For channel (2) the largest shift was +10%. 

To estimate the uncertainties in the tracking procedure, all track selection cuts were 
varied. In the present analysis only tracks associated with the event vertex were con- 
sidered. The systematic error due to any D* track not being fitted to the vertex was 
estimated from a careful comparison of vertex fitting in data and Monte Carlo. The re- 
sulting uncertainties on the cross sections are 1l\ 4 % and t\o% for channels (1) and (2), 
respectively. 

The uncertainty in the acceptance of the E e > 10 ° > 12 GeV cut was estimated by 
changing the cut value by ±0.5 GeV. A variation in the cross section of t % for channel (1) 
and ty% for channel (2) was found. A shift of ±3% in the calorimeter energy scale 
produces a variation of for channel (1) and t?% for channel (2). 

Background estimation uncertainties were determined by varying the normalization 
region. The uncertainties were found to be 1?% for both channels (1) and (2). 

The parton density parametrizations [g§ MRS(G), MRS(A'),GRV94 HO and CTEQ3M 
for the proton and GRV-G HO, GRV-G LO, LAC-G1, GS-G HO and DG-G1 for the pho- 
ton were used in the MC event simulation. No significant variation of the acceptance was 
found in the kinematic region used for this analysis. 

Finally, contributions of ±1.5% from the luminosity measurement, and ±3.7% or 
±5.7% from the branching ratios [10| of channels (1) or (2) respectively, were included. 
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All contributions to the systematic errors were added in quadrature. The final sys- 
tematic errors to the total cross sections are given in table 1. For the differential cross 
sections they were added in quadrature to the statistical errors and are indicated as the 
outer error bars in Figures 3 and 4. 



6 Comparison with NLO QCD Calculations 



Full NLO computations of differential cross sections for heavy quark production in any 
kinematic region |J[ [| became available recently. We compare our measurements with 
two such calculations. 

One of the two computations was done in the massive charm scheme |§, where m c 
acts as a cutoff for the perturbative calculation. The program for cc photoproduction [p9f 
from this computation was used to produce total and differential D* cross sections in 
the restricted kinematic region of our measurement. From MC studies it was found 
that the contribution of bb production to the D* cross section is below 5%. The charm 



fragmentation into D* was performed using the Peterson formula ||30|| : 



/(*)« 




(1-*), 



with e c = 0.06 [[H[] and using the branching ratio Br(c — > D ) = 0.260 ± 0.021 as 
measured by the OPAL collaboration ]32|. Here z is the fraction of the charm quark 
momentum taken by the D*. In order to convert the photoproduction differential cross 
sections into electroproduction cross sections, the W range from 115 to 280 GeV was 
divided into 15 GeV wide bins. The photoproduction cross sections were calculated in 
the center of each bin and the results were added after multiplying by the proper flux 
factors ||. 

A reference calculation was performed with the MRSG and GRV-G HO parton density 
parametrizations for the proton and photon, respectively. The renormalization scale used 
was /iR = m_i_ = \Jm1 + p]_ (m c = 1.5 GeV) and the factorization scales of the photon 
and proton structure functions was y,p — The predicted reference cross sections thus 
obtained are compared to the measured ones in Figures 3 and 4. Varying m c between 
1.2-1.8 GeV or ijlr between 0.5-2.0m_i_ changes the cross sections in our kinematic range 
within ±20%. Decreasing e c to 0.035 [32] increases the cross section by 15%. Using 
different proton or photon parton density parametrizations changes the cross section by 
less than 10%. All of the cross sections predicted by || using the reference parameters 
are lower than the data and the cross section in the kinematic region of the Kir channel 
is about 50% of the measured one. Better agreement with the data can be obtained with 
the choice of the parameters (jlr = 0.5 m± and m c = 1.2 GeV as shown in Figures 3 and 4. 
With the present statistics the shapes of the NLO predictions are in reasonable agreement 
with the data. 
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Recently, another type of NLO calculation || was compared to the HI and ZEUS 
preliminary results. In this approach the charm quark is treated according to the mass- 
less factorization scheme, which assumes charm to be one of the active flavours inside 
the proton and the photon, in contrast to the massive charm scheme. The differential 
distributions obtained with the massless approach NLO calculation are shown in Fig- 
ures 3 and 4. The parton density parametrizations used for the proton and photon were 
CTEQ4M []33 and GRV HO ^7j, respectively. The renormalization and factorization 
scales, as well as the values of m c and e c used are the same as in the reference calcu- 
lation of the massive-charm approach. The agreement with our data is good. Using a 
different proton parametrization (MRSG PEfl ) hardly changes the results while the pho- 
ton parametrization ACFGP-mc |34] reduces the cross section by 20% in our kinematic 
range |R5 . 



7 Summary 

The integrated and differential D* photoproduction cross sections in ep collisions at 
HERA have been measured with the ZEUS detector in the restricted kinematic region 
Q 2 < 4GeV 2 , 115 < W < 280 GeV, p?* > 3GeV and -1.5 < r] D * < 1.0. The cross 
section a ep ^D*x obtained using the channel D* + — > D°ng — > (K~n + )'Kg was measured 
to be (10.6 ± 1.7(stat.) ±J;jj (syst.))nb. Another D* decay channel, D* + -> D°n^ -> 
(K~7t + 7t + 7t~)7Ts , has been studied and good agreement with the Kir channel has been 
found in the region of overlap (p®* > 4 GeV). A NLO perturbative QCD massive charm 
scheme calculation predicts cross sections smaller than our measured values. Another 
NLO calculation in which the charm quark is treated according to the massless factoriza- 
tion scheme is in agreement with the data. The shapes of the differential cross sections 
da /dp®* , da/drj D * and da/dW are reasonably reproduced by both models. 
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Figure 1: AM and M(Kir) distributions for channel (1). In (a), the AM distributions 
for right charge combinations from the M(Ktt) signal region (full circles) and for wrong 
charge combinations from the region 1.68 < M(Ktt) < 2.04 GeV (dashed histogram) are 
shown. In (b), the M(Ktt) distributions for the right charge combinations from the AM 
signal region (full circles) and for wrong charge combinations from the AM signal region 
(dashed histogram) are shown. The full lines are the results of fits to a sum of a Gaussian 
and the functional form A ■ (AM - m^Y for (a) and exp(C + D ■ M(Kir)) for (b). 
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Figure 2: AM and M^Knnn) distributions for channel (2). In (a), the AM distributions 
are shown for the right charge combinations from the M{Kt\-kti) signal region (full circles), 
and for the background estimate using the wrong charge combinations (dashed histogram) 
as described in the text. In (b), the M^Kmnr) distributions are shown for the right charge 
combinations from the AM signal region (full circles) and for the background estimate 
using the wrong charge combinations (dashed histogram) as described in the text. The full 
lines are the results of fits to a sum of a Gaussian and the functional form A ■ (AM-m,) B 
for (a) and P 2 (x) = C + Dx + Ex 2 , where x = M{Ktxtxtx) for (b). 
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Figure 3: The differential cross section da ep -,o*x/ 'dp®* of D* photoproduction, 
Q 2 < 4 GeV 2 , in the kinematic region 115 < W < 280 GeV and -1.5 < r] D * < 1.0. 
The experimental points are drawn at the positions of the average values of an expo- 
nential fit in each bin. The inner part of the vertical error bars shows the statistical 
error, while the outer one shows the statistical and systematic errors added in quadra- 
ture. The HI points |7|] include statistical errors and systematical errors due to trigger 
efficiency The prediction of a NLO perturbative QCD calculation from a massive charm 
approach || is given by the dot-dashed curve, using MRSG and GRV-G HO as parton 
density parametrizations for the proton and photon respectively, fragmentation param- 
eter e c = 0.06, renormalization scale /ir = m± and m c = 1.5 GeV. The dotted curve is 
from the same calculation, but for = 0.5 and m c = 1.2 GeV. The full curve comes 
from the massless charm approach calculation ||, using the same parameters as for the 
dot-dashed curve, but with CTEQ4M taken for the parton density parametrization for 
the proton. 
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Figure 4: Differential cross sections for Q 2 < 4 GeV 2 : da ep ^o*x / di] D * in the kinematic 
region 115 < W < 280 GeV and pf* > 3 GeV (upper plot) and da ep -^d*x 1 1 dW in the 
kinematic region pj* > 3 GeV and —1.5 < r] D * < 1.0 (lower plot) for channel (1). The 
points are drawn at the centres of the corresponding bins. The inner part of the vertical 
error bars shows the statistical error, while the outer one shows the statistical and system- 
atic errors added in quadrature. The prediction of a NLO perturbative QCD calculation 
from a massive charm approach || is given by the dot-dashed curve, using MRSG and 
GRV-G HO as parton density parametrizations for the proton and photon respectively, 
fragmentation parameter e c = 0.06, renormalization scale hr = m± and m c = 1.5 GeV. 
The dotted curve is from the same calculation, but for fiR = 0.5 m± and m c = 1.2 GeV. 
The full curve comes from the massless charm approach calculation 0, using the same 
parameters as for the dot-dashed curve, but with CTEQ4M taken for the parton density 
parametrization for the proton. 
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